home *** CD-ROM | disk | FTP | other *** search
/ C/C++ Users Group Library 1996 July / C-C++ Users Group Library July 1996.iso / listings / v_02_12 / smith1.lis < prev    next >
Encoding:
File List  |  1991-09-24  |  346 b   |  14 lines
  1.  
  2. #include <math.h>
  3. #define ABSOLUTE_EPSILON      1e-10
  4. #define RELATIVE_EPSILON      1e-6   /* 6 digits */
  5. float xnew, xold;
  6.  
  7. xnew = 1.0;  /* some initial estimate */
  8. do { 
  9.   xold = xnew;
  10.   xnew = solution_refinement (xold);
  11. } while ( fabs(xnew - xold) > ABSOLUTE_EPSILON &&
  12.      fabs (xnew-xold) > fabs(xnew) * RELATIVE_EPSILON );
  13.  
  14.